Yield criteria for quasibrittle 
and frictional materials 



Davide Bigoni and Andrea Piccolroaz 
Dipartimento di Ingegneria Meccanica e 

Strutturale, Universita di Trento, 
Via Mesiano 77, 1-38050 Trento, Italia 
email: bigoni@ing.unitn.it 

August 29, 2002 

Abstract 

A new yield/ damage function is proposed for modelling the inelastic behaviour of 
a broad class of pressure-sensitive, frictional, ductile and brittle-cohesive materials. 
The yield function allows the possibility of describing a transition between the shape 
of a yield surface typical of a class of materials to that typical of another class of 
materals. This is a fundamental key to model the behaviour of materials which be- 
come cohesive during hardening (so that the shape of the yield surface evolves from 
that typical of a granular material to that typical of a dense material) , or which de- 
crease cohesion due to damage accumulation. The proposed yield function is shown 
to agree with a variety of experimental data relative to soil, concrete, rock, metallic 
and composite powders, metallic foams, porous metals, and polymers. The yield 
function represents a single, convex and smooth surface in stress space approaching 
as limit situations well-known criteria and the extreme limits of convexity in the 
deviatoric plane. The yield function is therefore a generalization of several crite- 
ria, including von Mises, Drucker-Prager, Tresca, modified Tresca, Coulomb-Mohr, 
modified Cam-clay, and — concerning the deviatoric section — Rankine and Ot- 
tosen. Convexity of the function is proved by developing two general propositions 
relating convexity of the yield surface to convexity of the corresponding function. 
These propositions are general and therefore may be employed to generate other 
convex yield functions. 

Keywords: B. yield criteria; B. elastic-plastic material; B. concrete; B. foam material; 
geological material; B. granular material 



1 Introduction 



Yielding or damage of quasibrittle and frictional materials (a collective denomination for 
soil, concrete, rock, granular media, coal, cast iron, ice, porous metals, metallic foams, 
as well as certain types of ceramic) is complicated by many effects, including dependence 
on the first and third stress invariants (the so-called 'pressure-sensitivity' and 'Lode- 
dependence' of yielding), and represents the subject of an intense research effort. Re- 
stricting attention to the formulation of yield criteria, research moved in two directions: 
one was to develop such criteria on the basis of micromechanics considerations, while 
another was to find direct interpolations to experimental data. Examples of yield func- 
tions generated within the former approach are numerous and, as a paradigmatic case, 
we may mention the celebrated Gurson criterion (Gurson, 1977). The latter approach 
was also broadly followed, providing some very successful yield conditions, such as for 
instance the Ottosen criterion for concrete (Ottosen, 1977). Although very fundamental 
in essence, the micromechanics approach has limits however, particularly when employed 
for geomaterials. For instance, it is usually based on variational formulations, possible 
— for inelastic materials — only for solids obeying the postulate of maximum dissipation 
at a microscale, which is typically violated for frictional materials such as for instance 
soils. 

A purely phenomenological point of view is assumed in the present article, wherein 
a new yield function^ is formulated, tailored to interpolate experimental results for qua- 
sibrittle and frictional materials, under the assumption of isotropy. The interest in this 
proposal lies in the features evidenced by the criterion. These are: 

• finite extent of elastic range both in tension and in compression; 

• non-circular deviatoric section of the yield surface, which may approach both the 
upper and lower convexity limits for extreme values of material parameters; 

• smoothness of the yield surface; 

• possibility of stretching the yield surface to extreme shapes and related capability 
of interpolating a broad class of experimental data for different materials; 

• reduction to known criteria in limit situations; 

• convexity of the yield function (and thus of the yield surface); 

• simple mathematical expression. 

1 We need not distinguish here between yield, damage and failure. Within a phenomenological 
approach, all these situations are based on the concept of stress range, bounded by a given hypersurface 
defined in stress space. 
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None of the above features is essential, in the sense that a plasticity theory can be 
developed without all of the above, but all are desirable for the development of certain 
models of interest, particularly in the field of geomaterials. This is a crucial point, 
deserving a carefully explanation. In particular, while some of the above requirements 
have a self-evident meaning, smoothness and convexity need some discussion. 

Although experiments are inconclusive in this respect (Naghdi et al. 1958; Paul, 
1968; Phillips, 1974), theoretical speculations (sometimes criticized, Naghdi and Srini- 
vasa, 1994) suggest that corners should be expected to form in the yield surface for single 
crystals and polycrystals (Hill, 1967). Therefore, smoothness of the yield surface might 
be considered a mere simplification in the constitutive modelling of metals. However, 
the situation of quasibrittle and frictional materials is completely different. For such 
materials, in fact, evidence supporting corner formation is weakcl, so that, presently, 
smoothness of the yield surface is a broadly employed concept and models developed 
under this assumption are still very promising. Moreover, corners often are included in 
the constitutive description of a material for the mere fact that an appropriate, smooth 
yield function is simply not available (this is usually the case of the apex of the Drucker- 
Prager yield surface and of the corner which may exist at the intersection of a smooth, 
open yield surface with a cap). 

Regarding convexity of the yield surface, we note that this follows for polycrystals 
from Schmid laws of single crystals (Bishop and Hill, 1951; Mandel, 1966). However, 
differently from smoothness, convexity is supported by experiments in practically all 
materials and is a useful mathematical property, which is the basis of limit analysis 
and becomes of fundamental importance in setting variational inequalities for plasticity 
(Duvaut and Lions, 1976). We may therefore conclude that — in the absence of a clear and 
specific motivation — it is not sensible to employ a yield function that violates convexity. 

A number of failure surfaces have been proposed meeting some of the above require- 
ments, among others, we quote the Willam and Warnke (1975), Ottosen (1977) and 
Hisieh et al. (1982) criteria for concrete, the Argyris et al. (1974), Matsuoka and Nakai 
(1977), Lade and Kim (1995), and Lade (1997) criteria for soils. For all these criteria, 
while some information can be found about the range of parameters corresponding to 
convexity of the yield surfaces, nothing is known about the convexity of the corresponding 
yield functions. 

Convexity of a yield function implies convexity of the corresponding yield surface, 
but convexity of a level set of a function does not imply convexity of the function itself. 
While it can be pointed out that a convex yield function can in principle always be found 
to represent a convex yield surface, the 'practical problem' of finding it in a reasonably 
simple form may be a formidable one. From this respect, general propositions would be 
of interest, but the only contribution of which the authors are aware in this respect is 
quite recent (Mollica and Srinivasa, 2002). A purpose of the present paper is to pro- 

2 Some argument in favour of corner formation in geomaterials have been given by Rudnicki and Rice 
(1975). 
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vide definitive results in this direction. In particular, the range of material parameters 
corresponding to convexity of the yield function proposed in this paper is obtained by 
developing a general proposition that can be useful for analyzing convexity of a broad 
class of yield functions. The proposition is finally extended to introduce the possibility of 
describing a modification in shape of the deviatoric section with pressure. The proposi- 
tions are shown to be constructive, in the sense that these may be employed to generate 
convex yield functions (examples of which are also included). 

Beyond the issue of convexity, the central purpose of this paper is the proposal of 
a yield criterion [see eqns. ©-©]• This meets all of the above-listed requirements 
and can be viewed as a generalization of the following criteria: von Mises, Drucker- 
Prager, Tresca, modified Tresca, Coulomb-Mohr, modified Cam-clay, Deshpande and 
Fleck (2000), Rankine, and Ottosen (1977) (the last two for the deviatoric section). 
Obviously, the criterion may account for situations which cannot be described by the 
simple criteria to which it reduces in particular cases. Several examples of this may be 
found in the field of granular media, where several ad hoc yield conditions have been 
proposed, which may describe one peculiar material, but cannot describe another. In the 
present paper, it is shown with several examples that our yield criterion provides a unified 
description for a extremely broad class of quasi-brittle and frictional materials. Beyond 
the evident interest in generalization, there is a specific motivation for advocating the 
necessity of having a single criterion describing different materials. This lies in the fact 
that during hardening, a yield surface may evolve from the shape typical of a certain 
material to that typical of another. An evident example of this behaviour can be found in 
the field of granular materials, referring in particular to metal powders. These powders 
become cohesive during compaction, so that the material is initially a true granular 
material, but becomes finally a porous metal, whose porosity may be almost completely 
eliminated through sintering. The key to simulate this process is plasticity theory, so that 
a yield function must be employed evolving from the typical shape of a granular material 
('triangular' deviatoric and 'drop-shaped' meridian sections), to that of a porous metal 
(circular deviatoric and elliptic meridian sections) and, in case of sintering, to that of a 
fully-dense metal (von Mises criterion). Another example of extreme shape variation of 
yield function during hardening is the process of decohesion of a rock-like material due to 
damage accumulation, a situation in a sense opposite to that described above. Evidently, 
a continuous distortion of the yield surface can be described employing the criterion 
proposed in this paper and simply making material parameters depend on hardening. 



2 A premise on Haigh-Westergaard representation 

The analysis will be restricted to isotropic behaviour, therefore the Haigh-Westergaard 
representation of the yield locus is employed (Hill, 1950). This is well-known, so that we 
limit the presentation here to a few remarks that may be useful in the following. First, 
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we recall that: 

Al. a single point in the Haigh-Westergaard space is representative of the infinite (to 
the power three) stress tensors having the same principal values; 

A2. due to the arbitrariness in the ordering of the eigenvalues of a tensor, six different 
points correspond in the Haigh-Westergaard representation to a given stress tensor. 
As a result, the yield surface results symmetric about the projections of the principal 
axes on the deviatoric plane (Fig. CP; 

.A3, the Haigh-Westergaard representation preserves the scalar product only between 
coaxial tensors; 

AA. a convex yield surface — for a material with a fixed yield strength under triax- 
ial compression — must be internal to the two limit situations shown in Fig. [1] 
(Haythornthwaite, 1985). Note that the inner bound will be referred as 'the Rank- 
ine limit'. 

Due to isotropy, the analysis of yielding can be pursued fixing once and for all a reference 
system and restricting to all stress tensors diagonal in this system. We will refer to this 
setting as to the Haigh-Westergaard representation. When tensors (for instance, the yield 
function gradient) coaxial to the reference system are represented, the scalar product is 
preserved, property A3. In the Haigh-Westergaard representation, the hydrostatic and 
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Figure 1: Deviatoric section: definition of angle 9, symmetries, lower and upper convexity bounds. 



deviatoric stress components are defined by the invariants 



where 

J 2 = Vs, S = <x-^I, (2) 

in which S is the deviatoric stress, I is the identity tensor, a dot denotes scalar product 
and tr denotes the trace operator, so that A ■ B = tr AB T , for every second-order tensors 
A and B. The position of the stress point in the deviatoric plane is singled out by the 
Lode (1926) angle 9 defined as 

= ± cos" 1 (^^rV J 3 = -trS 3 , (3) 
3 ^2 J 2 3 / 2 j ' 6 3 V ; 

so that 9 G [0,7r/3]. As a consequence of property (.4.2) of the Haigh-Westergaard 
representation, a single value of 9 corresponds to six different points in the deviatoric 
plane (Fig. [TJ. The following gradients of the invariants, that will be useful later, 

*=_Ii ^i = s dJs -s 2 trS \ 

da 3 ' da ' da 3 

(4) 

m_ 9 / 2 trS 2 cos3^ s N 

da~ ~ 2 q 3 sin 36 \ 3~ ~ 9 3 

can be obtained from well-known formulae (e.g. Truesdell and Noll, 1965, Sect. 9) using 
the identity 

<9S 1 

_=I®I--I®I, (5) 
oa — 3 

where the symbol ® denotes the usual dyadic product and I® I is the symmetrizing 
fourth-order tensor, defined for every tensor A as I®I[A] = (A + A T )/2. Note that 
89 /da is orthogonal to I and to the deviatoric stress S. 



3 A new yield function 

We propose the seven-parameters yield function F : Sym — > K U {+00} defined as: 

F(<r) = m + WY (6) 

where the dependence on the stress a is included in the invariants p, q and 9, eqns (pQ) 
and ([3]), through the 'meridian' function 



-Mp cV /($ - $ m ) [2(1 - a)$ + a] if $ G [0, 1], 

f{p) = { (7) 
+00 i/$^[0,l], 
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where 

* = 4^, (8) 

p c + C 

describing the pressure-sensitivitjH and the 'deviatoric' function 

9(8) = 1 -ZI7 57777' ( 9 ) 

cos [p~ — I cos 1 (7 cos 3#)J 

describing the Lode-dependence of yielding. The seven, non-negative material parame- 
ters: 

M > 0, p c > 0, c > 0, < a < 2, m > 1, < /3 < 2, < 7 < 1, (10) 

V v ' S v ' 

defining f(p) defining g{9) 

define the shape of the associated (single, smooth) yield surface. In particular, M controls 
the pressure-sensitivity, p c and c are the yield strengths under isotropic compression and 
tension, respectively. Parameters a and m define the distortion of the meridian section, 
whereas ft and 7 model the shape of the deviatoric section. Note that the deviatoric 
function describes a piecewise linear deviatoric surface in the limit 7 — > 1. Finally, it 
is important to remark that within the interval of ft e [0, 2] the yield function is convex 
independently of the values assumed by parameter 7. Convexity requirements, that will 
be proved later, impose a broader variation of ft than (110[) k. but the interval where ft 
may range becomes a function of 7. In particular, the yield function is convex when 

2-B( T )</3<B( 7 ), (11) 

where function -6(7) takes values within the interval ]2,4], when 7 ranges in [0, 1[ and is 
defined as 

6 , 1 - 2 cos 2 - 2 cos 2 2: 
0(7) = 3 tan 



(12) 

2 = 2/3(7r — cos -1 7) 



7r 2sinz(l — cosz) 

The yield function ((6]) corresponds to the following yield surface: 

q = -f(p)g(0), pe[-c,p c ], ee[o,n/3], (13) 

which makes explicit the fact that f(p) and g(9) define the shape of the meridian and 
deviatoric sections, respectively. 

The yield surface f|T3|) is sketched in Figs. |2]|3] for different values of the seven above- 
defined material parameters (non-dimensionalization is introduced through division by 
p c in Fig. |2]). In particular, meridian sections are reported in Figs. [2] (g(9) = 1 has been 
taken), whereas Fig. [3] pertains to deviatoric sections. 



3 Thc meridian function can be written in an alternative form by using the Macauley bracket operator, 
denned for every scalar a as < a >= max{0, a}, and the indicator function X[o.i]( < ^): which takes the 
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Figure 2: Meridian section: effects related to the variation of parameters M (a), c/p c (b), m (c), and 
a(d). 



As a reference, the case corresponding to the modified Cam-clay introduced by Roscoe 
and Burland (1968) and Schofield and Wroth (1968) and corresponding to /3 — 1, 7 = 0, 
a — 1, m = 2, and c = is reported in Fig. [2] as a solid line, for M = 0.75. The 
distortion of meridian section reported in Fig. [2] (a) — where M = 0.25, 0.75, 1.25 — can 
also be obtained within the framework of the modified Cam-clay, whereas the effect of an 
increase in cohesion reported in Fig. [2] (b) — where c/p c = 0, 0.2, 0.4 — may be employed 
to model the gain in cohesion consequent to plastic strain, during compaction of powders. 

The shape distortion induced by the variation of parameters m and a, Fig. [2] (c) - 
where m = 1.2, 2, 4 — and (d) — where a = 0.01, 1.00, 1.99 — is crucial to fit experimental 
results relative to frictional materials. 

A unique feature of the proposed model is the possibility of extreme shape distortion 
of the deviatoric section, which may range between the upper and lower convexity limits, 
and approach Tresca, von Mises and Coulomb-Mohr. This is sketched in Fig. [31 where 
to simplify reading of the figure, function g{9) has been normalized through division 
by g(ir/3), so that all deviatoric sections coincide at the point 9 = ir/3. The use of 



value when $ 6 [0, 1] and is equal to +00 otherwise 

f(p) = -Mp c J($-$ m ) 2(1 - a)$ + a + X[o,i]($), $ =<$>-<$-!> 



(a) y = 0.99 s (b) y = 1 




<c) P = 9 {d} p = 0.5 9 




Figure 3: Deviatoric section: effects related to the variation of (3 and 7. Variation of ft = 0, 0.5, 1, 1.5, 2 
at fixed 7 = 0.99 (a) and 7 = 1 (b). Variation of 7 = 1, 0.75, at fixed /3 = (c) and (3 = 0.5 (d). 



our model may therefore allow one to simply obtain a convex, smooth approximation of 
several yielding criteria (Tresca and Coulomb- Mohr, for instance). If this may be not 
substantial from theoretical point of view, it clearly avoids the necessity of introducing 
independent yielding mechanisms. 

Parameter 7 is kept fixed in Figs. |3] (a) and (b) and equal to 0.99 and 1, respectively, 
whereas parameter (3 is fixed in Figs. [3] (c) and (d) and equal to and 1/2. Therefore, 
figures (a) and (b) demonstrate the effect of the variation in /3 (= 0, 0.5, 1, 1.5, 2) which 
makes possible a distortion of the yield surface from the upper to lower convexity limits 
going through Tresca and Coulomb-Mohr shapes. The role played by 7 (= 1, 0.75, 0) is 
investigated in figures (c) and (d), from which it becomes evident that 7 has a smoothing 
effect on the corners, emerging in the limit 7 = 1. The von Mises (circular) deviatoric 
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section emerges when 7 = 0. 

The yield surface in the biaxial plane o\ versus a 2 , with <t 3 = is sketched in Fig. HJ 
where axes are normalized through division by the uniaxial tensile strength f t . In par- 
ticular, the figure pertains to M = 0.75, p c = 50c, m = 2, and a = 1, whereas 7 = 0.99 
is fixed and (3 is equal to {0, 0.5, 1, 1.5, 2} in Fig. @] (a) and, vice-versa, (3 = is fixed 
and 7 is equal to {0, 0.75, 0.99} in Fig. H (b). 



(a) y = 0.99 



<b) p = 



,11 



a, /f t 




tr, II 




a, ft 



Figure 4: Yield surface in the biaxial plane tJ\j ft vs. (T2/ ft, with 03 = 0. Variation of / 
at fixed 7 = 0.99 (a) and variation of 7 = 0, 0.75, 0.99 at fixed /3 = (b). 



0,0.5,1,1.5,2 



3.1 Smoothness of the yield surface 

Smoothness of yield surface (fT3|) within the interval of material parameters defined in 
( fl~0l) - (fTTl) can be proved considering the yield function gradient. This can be obtained 
from (jll) in the form 

BF 

a(p)l + b(9)S + c(9)S L , (14) 



No- 



where 

and 
a{p) 

b(9) 
c{9) 



3S - ± _ V2 09 
2? ' 

19/(p) 

3 dp 3(p c + c) 

3 1 



cos 3$ S 



y/Sqder sin 39 

Mp c (1 - mr- 1 ) [2(1 - a)$ + a] + 2(1 - a) ($ - $ r 



(15) 



2 v /($ _ $m) [2(1 - a )$ + 



2^)' 

\/37 sin 36* 
V2^l - 7 2 cos 2 30 



(16) 



sm 



7T 1 

/3 cos -1 (7 cos 39) 

6 3 
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It should be noted that c(0) = c(n/3) = and that S and S -1 are unit norm, coaxial and 
normal to each other tensors^. Coaxiality and orthogonality are immediate properties, 
whereas the proof that | !S — | = 1 is facilitated when the following identities are kept into 
account 

~ 3 _ 1 - _ cos30 ^ s 2 -S 2 = i, (17) 

the former of which is the Cayley-Hamilton theorem written for S. Let us consider now 
from (j!4)) the unit-norm yield function gradient 

Q= ; a =I+- S + ^= ° S x , (18) 

VSa 2 + b 2 + c 2 V3a 2 + b 2 + c 2 V3a 2 + b 2 + c 2 

defining, for stress states satisfying F(cr) = 0, the unit normal to the yield surface. The 
following limits can be easily calculated 

lim Q = -i=I, lim Q = — ^1, (19) 

so that the yield surface results to be smooth at the limit points where the hydrostatic 
axis is met. Moreover, smoothness of the deviatoric section of the yield surface is proved 
observing that 

lim Q= - a -A + - b -_S, (20) 
e-*o,ir/3^ V3a 2 + b 2 V3a 2 + b 2 

where S and b are evaluated at 9 = and 9 = tt/3, and noting that S and S -1 are 
coaxial, deviatoric tensors so that they are represented by two orthogonal vectors in the 
deviatoric plane in the Haigh-Westergaard stress space. We observe, finally, that limits 
(fl9~|) do not hold true when a equals and 2 and that limits ( 1201) does not hold true when 
7 = 1. In particular, a corner appears at the intersection of the yield surface with the 
hydrostatic axis in the former case and the deviatoric section becomes piecewise linear 
in the latter. 



3.2 Reduction of yield criterion to known cases 

The yield function (JHJ)-(IH1) reduces to almost aljfl 'classical' criteria of yielding. These can 
be obtained as limit cases in the way illustrated in Tab. [1] where the modified Tresca 
criterion was introduced by Drucker (1953), whereas the Haigh-Westergaard representa- 
tion of the Coulomb-Mohr criterion was proposed by Shield (1955). In Tab. [1] parameter 

4 Note that cS 1 - = at 6 = 0,ir/3. This can be deduced from the fact that | S ~ *- 1 = 1 and c = for 
9 = 0,7r/3 or, alternatively, can be proved directly observing that for 9 = 0,7r/3 the deviatoric stress 
can be generically written as {Si, —S±/2, —5*1/2}, unless all (uninfluent) permutations of components. 

5 A remarkable exception is the isotropic Hill (1950 b) criterion, corresponding to a Tresca criterion 
rotated of ir/6 in the deviatoric plane. 
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r denotes the ratio between the uniaxial strengths in compression (taken positive) and 
tension, indicated by f c and f t , respectively. We note that for real materials r > 1 and 
that we did not explicitly consider the special cases of no-tension f t = or granular 
ft = fc = materials [which anyway can be easily incorporated as limits of 

We note that the expression of the Tresca criterion which follows from ([6])- (Q in the 
limits specified in Tab. [I] was provided also by Bardet (1990) and answers — in a positive 
way — the question (raised by Salengon, 1974) if a propel form of the criterion in terms 
of stress invariants exist. 

The Mohr-Coulomb limit merits a special mention. In fact, if the following values of 
the parameters are selected 

f c [cos (M - f ) + cos Ml 3 \r cos (M - \ ) - cos Ml , 

a = 0, c=-^ ^ V — , M = -± K —l= — ^ 21 

3rcos(/3f -f) -3cos/3f y/2{r + 1) 

and then the limits 

7 — ► 1, Pc = f c m — y oo, (22) 

are performed, a three-parameters generalization of Coulomb-Mohr criterion is obtained, 
which reduces to the latter criterion in the special case when /3 is selected in the form 
specified in Tab. [T] (yielding an expression noted also by Chen and Saleeb, 1982). The 
cases reported in Tab. prefer to situations in which the criterion (EJ)- ([9j) reduces to known 
yield criteria both in terms of function f(p) and of function g(6). It is however important 
to mention that the Lode's dependence function g{9) reduces also to well-known cases, 
but in which the pressure-sensitivity cannot be described by the meridian function (J?]). 
These are reported in Tab. [2j It is important to mention that the form of our function 
g(9), eqn. flS]), was indeed constructed as a generalization of the deviatoric function 
introduced by Ottosen (1977). 



3.3 A comparison with experiments 

A brief comparison with experimental results referred to several materials is reported 
below. We limit the presentation to a few representative examples demonstrating the 
extreme flexibility of the proposed model to fit experimental results. In particular, we 
concentrate on the meridian section, whereas only few examples are provided for the 
deviatoric section, which has a shape so deformable and ranging between well-known 

6 The expression 

/(<t) = 4J| - 27 J 3 2 - 36k 2 J\ + 96fc 4 J 2 - 64/c 6 , 

where k is the yield stress under shear (i.e. k = ft/2), reported in several textbooks on plasticity, is 
definitively wrong. This can be easily verified taking a stress state belonging to one of the planes defining 
the Tresca criterion, but outside the yield locus, for instance, the point {a\ =0, 02 = — 2fc, 03 = 2k}, 
corresponding to J2 = 4fc 2 and J3 = 0. Obviously, the point lies well outside the yield locus, but satisfies 
/(er) = 0, when the above, wrong, yield function is used. 
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Table 1: Yield criteria obtained as special cases of ([(J])-©, r — fd ft an d f c and f t are the uniaxial 
strengths in compression and tension, respectively. 



Criterion 


Meridian function f(p) 


Deviatoric function 


von Mises 


a — 1 m = 2, 
i\/r ^ft 

M = — , c = p c = — > oo 

Pc 


£ = 1,7 = 


Drucker-Prager 


a = 0, M — —= , 

V2(r + 1) 

c = o ( 1V Pc = f c m — > oo 
3(r — 1) 


as for von Mises 


Tresca 


as for von Mises, except that 
M=^ ft 

Pc 


= 1,7— >1 


mod. Tresca 


as for Drucker-Prager, except that 
M= 3V3(r-U 


as for Tresca 


Coulomb-Mohr 


as for Drucker-Prager, except that 
M 3[rcos(/3f-f)-cos/3|] 


a 6 j V3 

p = — tan , 


V^r + 1) 
/ c [cos(/3|-f)+cos/9|] 
3rcos(/3f - f) -3cos/3f 


7T 2r + 1 ' 


mod. Cam-clay 


m = 2, a = 1, c = 


as for von Mises 
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Table 2: Deviatoric yield functions obtained as special cases of © 



Criterion 


Deviatoric function g{6) 


Lower convexity (Rankine) 


P = 0, 7— ►! 


Upper convexity 


(3 = 2, T — U 


Ottosen 


(3 = 0, < 7 < 1 



forms that fitting experiments is a-priori expected. Results on the biaxial plane a% — a 2 
are also included. All values of material parameters defining the yield function ([S])-© 
employed to fit experimental data may be useful as a reference and are reported in 
Appendix A. 

Typical of soils are the experimental results reported in Fig. |5l on Aio dry sand and 
Weald clay, taken, respectively, from Yasufuku et al. (1991, their Fig. 10a) and Parry 
(reported by Wood, 1990, their Fig. 7.22, so that p e is the equivalent consolidation pres- 
sure in Fig. E£b)). Note that the upper plane of the graphs refers to triaxial compression 
(6 = 7r/3), whereas triaxial extension is reported in the lower part of the graphs (6 = 0). 
It may be concluded from the figure that experimental results can be easily fitted by our 
function f(p), still maintaining a smooth intersection of the yield surface with p— axis. 

(a) (b) 

°' 8 Q = n/3 6 = 7i/3 




Figure 5: Comparison with experimental results relative to sand (Yasufuku et al. 1991) (a) and clay 
(Parry, reported by Wood, 1990) (b). 



In addition to soils, the proposed function can model yielding of porous ductile 

or cellular materials, metallic and composite powders, concrete and rocks. To further 
develop this point, a comparison with experimental results given by Sridhar and Fleck 
(2000) — their Figs. 5(b) and 9(c) — relative to ductile powders is reported in Fig. [6j In 
particular, Fig. |6](a) is relative to an aluminum powder (Al Do=0.67, D=0.81 in Sridhar 
and Fleck, their Fig. 5b), Fig. |6] (b) to an aluminum powder reinforced by 40 vol.%SiC 
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(Al-40%Sic D =0.66, D=0.82 in Sridhar and Fleck, their Fig. 5b), Fig. El (c) to a lead 
powder (0% steel in Sridhar and Fleck, their Fig. 9c), and Fig. El(d) to a lead shot-steel 
composite powder (20% steel in Sridhar and Fleck, their Fig. 9c). Beside the fairly good 
agreement between experiments and proposed yield function, we note that the aluminium 
powder has a behaviour — different from soils and lead-based powders — resulting in a 
meridian section of the yield surface similar to the early version of the Cam-clay model 
(Roscoe and Schofield, 1963). 



9 = 7t/3 




(b) 



e = ti/3 




0.2 0.4 0.6 0. 




Figure 6: Comparison with experimental results relative to aluminium powder (a) aluminum composite 
powder (b), lead powder (c) and lead shot-steel composite powder(d), data taken from Sridhar and Fleck 
(2000). 



Regarding concrete, among the many experimental results currently available, we 
have referred to Sfer et al. (2002, their Fig. 6) and to the Newman and Newman (1971) 
empirical relationship 

/ \ 0.86 

where 0\ and cr 3 are the maximum and minimum principal stresses at failure and f c is 
the value of the ultimate uniaxial compressive strength. Small circles in Fig. [7] represents 
results obtained using relationship (|2"3"|) in figure (a) and experimental results by Sfer 
et al. (2002) in figure (b); the approximation provided by the criterion ©-((8]) is also 
reported continuous line. 

As far as rocks are concerned, we limit to a few examples. However, we believe that 
due to the fact that our criterion approaches Coulomb-Mohr, it should be particularly 
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suited for these materials. In particular, data taken from Hoek and Brown (1980, their 
pages. 143 and 144) are reported in Fig. El as small circles for two rocks, chert (Fig. [8] a) 
and dolomite (Fig. El b). 



(a) 6 = 7i/3 (b) 6 = 7i/3 




0.511.52 0.5 11.52 



P / f c P / f c 

Figure 8: Comparison with experiments for rocks. Chert (a) and dolomite (b), data taken from Hoek 
and Brown (1980). 

A few data on polymers are reported in Fig. |9] — together with the fitting provided by 
our model — concerning polymethil methacrylate (Fig. [9] a) and an epoxy binder (Fig. [9] 
b), taken from Ol'khovik (1983, their Fig. 5), see also Altenbach and Tushtev (2001, 
their Figs. 2 and 3). 
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(a) 



B = n/3 



(b) 



Q = n/3 



CD 
CL 




100 200 300 400 

p (MPa) 




100 200 300 400 

p (MPa) 



Figure 9: Comparison with experimental results for polymers. Methacrylate (a) and an epoxy binder 
(b), data taken from Ol'khovik (1983). 



Finally, our model describes — with a different yield function — the same yield surface 
proposed by Deshpande and Fleck (2000) to describe the behaviour of metallic foams. In 
particular, the correspondence between parameters of our model (JSJ)- (J9j) and of the yield 
surface proposed by Deshpande and Fleck [2000, their eqns. (2)-(3)] is obtained setting 

{3 = 1, 7 = 0, m = 2, a = 1, 

and assuming the correlations given in Tab. [3j 



Table 3: Correspondence between parameters of (|6 ]) -(|9" ]l and Deshpande and Fleck (2002) yield functions 
— the latter shortened as 'DF model' — to describe the behaviour of metallic foams. 





Model Q6 




DF model 


M 


c 


Pc 
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DF model 

Y,a 


2a 










Model (JED-© 

M,p c ,c 








cM 


M 

T 








2^+(f) 2 



The proposed function (jSJ)-© is also expected to model correctly yielding of porous 
ductile metals. As a demonstration of this, we present in Fig. [10] a comparison with the 
Gurson (1977) model. The Gurson yield function has a circular deviatoric section so that 
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/3 — 1 and 7 = in our model, in addition, we select 
1, m = 2, p, 



a 



c = a M -^- cosh- 1 l+ J* f \ M = a M - v / l + g 3 / 2 -2/g 1 , (24) 

where / is the void volume fraction (taking the values {0.01, 0.1, 0.3, 0.6 } in Fig. [10]), &m 
is the equivalent flow stress in the matrix material and q% = 1.5, q<± = 1 and ^3 = q\ are 
the parameters introduced by Tvergaard (1981, 1982). A good agreement between the 
two models can be appreciated from Fig. UHl increasing when the void volume fraction / 
increases. 



1 

.8 
. 6 
0,4 
0,2 



Present model 













\q.i 




\0.3 




\ 1 
\ 1 


,0.6 \ 







0. 5 



1.5 2 



2. 5 



3. 5 



Figure 10: Comparison with the Gurson model at different values of void volume fraction f. 



As far as the deviatoric section is regarded, we limit to two examples — reported in 
Fig. [H] — concerning sandstone and dense sand, where the experimental data have been 
taken from Lade (1997, their Figs. 2 and 9a). 




30 SB 




Figure 11: Comparison with experimental results relative to deviatoric section for sandstone (a) and 
dense sand (b) data taken from (Lade, 1997). 



Experimental data referred to the biaxial plane a 3 = for grey cast iron and concrete 
(taken respectively from Coffin and Schenectady, 1950, their Fig. 5 and Tasuji et al. 
1978, their Figs. 1 and 2) are reported in Fig. [12j 
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(a) Bi/f B (b) cr 5 /f c 



0.5, 0.5 




-1.5 -1.5 



Figure 12: Comparison with experimental results on biaxial plane for cast iron (data taken from Coffin 
and Schenectady, 1950) (a) and concrete (data taken from Tasuji et al. 1978) (b). 



4 On convexity of yield function and yield surface 

Convexity of the yield function (jS])-® within the range of parameters (TTU|) - (ITT|) was 
simply stated in the previous Section and still needs a proof. This will be given at the 
end of the present section as an application of a general proposition relating convexity of 
yield functions and surfaces that is given below. 

We begin noting that while convexity of yield function implies convexity of the corre- 
sponding yield surface, the converse is usually false, namely, convexity of the level set of 
a function is unrelated to convexity of the function itself. As an example, let us consider 
the non- convex yield function 

f(p,q) = -,-~ 2 + q ^ 0<^<1, (25) 
a 4 a z b z a 

(where a and b are material parameters having the dimension of stress) which corresponds 
to a convex yield surface f(p,q) = 0, Fig. [131 After the pioneering work of de Finetti 
(1949), it became clear that convexity of every level set of a function represents its quasi- 
convexity, a property defining a class of functions much broader than the class of convex 
functions. In more detail, let us consider a function /(x) : U — > R, with U being a 
convex set, and its level sets 

L a = {xeU\ /(x) < a}, (26) 

so that: 

/ is quasi-convex if the level sets L a are convex for every a G R. 

Now, the above definition of quasi-convexity is equivalent (Roberts and Varberg, 1973) 
to the definition 

/[Ax+(l + A)y]<max{/(x),/(y)}, Vx,yGf/, VAg[0,1], (27) 
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and, if / is continuous and differentiable, to 

/(y)</(x) V/(x)-(x-y)>0, VxjGK (28) 

Convexity of the yield surface can be either accepted on the basis of experimental 
results, or on some engineering argumentation, such as for instance Drucker's postulate. 
Obviously, a convex yield locus can be expressed as a level set of a function, that generally 
may lack convexity and, even, quasi- convexity. For example, the level sets of function 
(|25|) are given in Fig. [13j It can be observed that while f(p, q) = may perfectly serve 




0.2 0.4 0.S 0.8 1 



p / a 

Figure 13: Level sets of function (|25|) . 



as a (convex) yield surface, the corresponding yield function even lacks quasi- convexit}d 
It is true that, in principle, a convex yield function can always be found to represent a 
convex yield surface, but to find this in a reasonably simple form may be an hard task. 
In other words, a number of yield functions that were formulated as an interpolation of 
experimental results still need a proof of convexity, even in cases where the corresponding 
yield locus is convex. The propositions that will be given below set some basis to provide 
these proofs. 

4.1 A general result for a class of yield functions 

The yield function (E])-® presented in the previous section may be viewed as an element 
of a family of models specified by the generic form (jH]). This family includes, among 
others, the models by Gudheus (1973), Argyris et al. (1974), Willam and Warnke (1975), 

7 As noted by Franchi et al. (1990), definition d28j is very similar to Drucker's postulate. However, 
Drucker's postulate merely prescribes the so-called normality rule of plastic flow and convexity of yield 
surface (Drucker, 1956, 1964). Quasi-convexity becomes a consequence of Drucker's postulate only in the 
special case — considered by Franchi et al. (1990) — in which convexity of yield surface implies convexity 
of all level sets of the corresponding function. 
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Eekelen (1980), Lin and Bazant (1986), Bardet (1990), Ehlers (1995), and Menetrey and 
Willam (1995) and Christensen (1997) and Christensen et al. (2002). 

A general result is provided below showing that for the range of material parameters 
for which the Haigh-Westergaard representation of a yield surface (fl3|) is convex, the 
function is also convex. 



Proposition 1: Convexity of the yield function ([6]) is equivalent to convexity of the 
meridian and deviatoric sections of the corresponding yield surface (fl3|) in the Haigh- 
Westergaard representation. In symbols: 



convexity of F(cr) = f{p) + 
where g(8) is a positive function. 



9(0) 



/">0, & 



9 



2g' 2 -gg">0, 



(29) 



Proof. It is a well-known theorem of convex analysis (Ekeland and Temam, 1976) that 
the sum of two convex functions is also a convex function. Since q, p and 9 are independent 
parameters, failure of convexity of f(p) or q/g(6) implies failure of convexity of F(er) and 
therefore convexity of both f(p) and q/g(6) are necessary and sufficient conditions for 
convexity of F(cr). 

Now, let us first analyze f{p)- The fact that convexity of f(p) as a function of a is 
equivalent to convexity of the meridian section follows from linearity of the trace operator, 
in view of the fact that p = — trcr/3. 

Second, the fact that convexity of qj ' g{&) as function of a is equivalent to the convexity 
of the deviatoric section follows from the 3 lemmas listed below. 

□ 

Lemma 1 (Hill, 1968): Convexity of an isotropic function of a symmetric (stress) tensor 
cr is equivalent to convexity of the corresponding function of the principal (stress) values 
(Tj (i = 1, 2, 3). In symbols, given: 



(01, 2 , 0- 3 ) 



then: 



A— -• Act > 
oa 



i=i 



(30) 



(31) 



where A denotes an ordered difference in the variables, so that, denoting with A and B 
two points in the tensor space 



A-^-Acr 
ocr 



fd4> 






o-a da- 



A-— A0j 

do; 




d(j) 




°i - °i ) , 



1,2,3. 



(32) 



(33) 
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Proof. That convexity of <fi implies convexity of </> is self-evident. The converse is not 
trivial and a proof was given by Hill (1968) with reference to an elastic strain energy 
function. The proof, omitted here for brevity, was later obtained also by Yang (1980) 
with explicit reference to a yield function. 

□ 

Lemma 2: Given a generic isotropic function of the stress that can be expressed as 

0(ai, 0- 2 , <t 3 ) = 4>(Si, S 2 ), (34) 
where Si and S 2 are two of the principal components of deviatoric stress, i.e. 

Si = - (2ai - a 2 - 0-3) , S 2 = - (-0-1 + 2cr 2 - cr 3 ) , (35) 

convexity of <p(ai, a 2 , cr 3 ) is equivalent to convexity of 4>(Si, S 2 )- 

Proof. The proof follows immediately from the observation that the relation ( 135]) 
between {Si, 62} and {o~i, a 2 , 0-3} is linear. 

□ 

Lemma 3: Convexity of 

m m 

as a function of Si, S 2 is equivalent to the convexity of the deviatoric section in the 
Haigh-Westergaard space: 

g 2 + 2g' 2 -gg">0. (37) 



Proof. The Hessian of ( |36l) is 



d 2 q/g(9) 
dStdSj 



-99' 



9' 



d 2 q 



dSidSj 



+ q(2g' 2 -gg" 



dq 89 dq 89 

+ ~7777~ 7777" + 9 



89 89 
] dSidSj 
8 2 9 



dSi 8Sj 



8Sj 8Si 



dSidSj 



(3? 



where % and j range between 1 and 2 and all functions q and 9 are to be understood as 
functions of Si and S 2 only. Derivatives of q may be easily calculated to be 



dq 



2Si - (-l)Vni, 



d 2 q 



27 



dS^Sj Aq 



(39) 
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where indices are not summed and vector m has the components 

{m} = {S 2 , -Si}. 

The derivatives of 9 can be performed through cos 3$, eqn. (j3])i, noting that 
89 



— 1 d cos 39 
8S~i = 3 sin 39 OSi ' 

d 2 9 -1 



cos 39 d cos 39 d cos 39 d 2 cos 36* 

+ 



so that 



dSidSj 3 sin 39 \ sin 2 39 OSi dSj dSidSjJ' 



dq 89 dq 89 d 2 9 



dSidSj dSjdSi * dS l dS j 
-1 



sin 39 



d 2 qcos39 d 2 q cos 39 d cos 39 d cos 39 

cos 39 - - — + g— 



sin 2 3^ 95, d^- 



where 



d cos 3# 9^ sin 39 d 2 qcos39 



rrii 



dSi 2q 2 
A substitution of <$S§ into Q32J) yields 

<9g <9# <9g <96> 



-27 — m^m,-. 
g 6 



<9 2 # 



asias* dftdft ^SidS 



so that we may conclude that the Hessian ( 1381) can be written as 



d 2 q/g(9) 27 (g 2 + 2g' 2 - gg"] 



dSidSj 



q 3 g 3 



rrii rrij . 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



Positive semi-definiteness of the Hessian (1451) is condition (1371) . which, in turn, represents 
non-negativeness of the curvature (and thus convexity) of deviatoric section. 

□ 



4.2 Applications of Proposition 1 

The scope of this section is on one hand to prove the convexity of function ([6])- (J9j) within 
the range ( !T0|) - (|TT|) of material parameters, on the other hand to show that Proposition 
1 is constructive, in the sense that can be used to invent convex yield functions. Let us 
begin with the first issue. 



23 



4.2.1 The proposed yield function 

First, we show that f(p), eqn (jBJ), is a convex function of p (so that the meridian section 
is convex) and, second, that the deviatioric section described by g(9), eqn (Q, is convex 
for the range of material parameters listed in (|T0l) - (|TT|) . Therefore, as a conclusion from 
Proposition 1, function F(er) results to be convex. 

A well-known result of convex analysis (Ekelan and Temam, 1976) states that function 
f(p) is convex if and only if the restriction to its effective domain (i.e. $ G [0,1]) is 
convex. Moreover, the function $ appearing in is a linear function of p so that 
convexity of f(p) can be inferred from convexity of the corresponding function, say /, of 
$. Introducing for simplicity the function 

= ($ - $ m ) [2(1 -a)$ + a], (46) 

the convexity of function /(<&) reduces to the condition 

[/i / ($)] 2 -2/i / '($)/i($) > 0, (47) 

where 

h'($) = (1 - m<$> m - 1 ) [2(1 - a)$ + a] + 2(1 - a) ($ - $ m ) , 

= -m(m - l)$ m - 2 [2(1 - a)$ + a] + 4(1 - a) (1 - m^™- 1 ) . 

Fulfillment of eqn f|4T|) can be now easily proven considering the inequality 

h"{$) <4(l-a) (1 - ml?" 1 - 1 ) , V$g[0,1]. (49) 

It remains now to show convexity of q/g(6). To this purpose, Proposition 1 can be 
employed, through substitution of © into the convexity condition eqn fl371) . thus yielding 



(48) 



1 37 cos 36* 

H . sin 



9(0) y/i - 7 2 cos 2 39 



IX 1 

- o cos_ (7 cos 30) 
6 3 



> 0, (50) 



where G [0, 7r/3] and g>(0) is given by eqn (Q. For values of 7 belonging to the interval 
specified in (TTUl) ^. condition (1501) can be transformed into 

sin - V) + 2 sin ^ + |^ > 0, (51) 

with x G [cos -1 7, 7r — cos -1 7] and then into 

— 1 + 2 COS 2; + 2 COS 2 Z „7T n , , 

sin/3- + COS/3- > 0, (52) 



2 sinz(l — cosz) 6 6 

with 2 G [2/3 cos -1 7, 2/3(7r — cos -1 7)], an inequality that can be shown to be verified 
within the interval of specified in (1TTT) and thus also within its subinterval (Pile- 
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4.2.2 Generating convex yield functions 



Proposition 1 can be easily employed to build convex yield functions within the class 
described by eqn <Q. The simplest possibility is to maintain f(p) in the form ([7]) and 
change the deviatoric function eqn (Q. As a first proposal, we can introduce the following 
function 

£(fl) = [l + /3(l + cos3#)r 1/n , (53) 

instead of . This describes a smooth deviatoric section approaching (without reaching) 
the triangular (Rankine) shape when parameters n > and > are varied. The yield 
function is convex within the range of parameters reported in Tab. H] (see Appendix B for 
a proof). The yield function defined by eqns ([7]) and (153]) does not possess the extreme 

Table 4: Conditions for the convexity of deviatoric yield function (|53|) . 



< n < 11/3 



n 



9-2n 



n > 11/3 



0< 



-1 + W1 + 



9(n - 2) 2 
n 2 (An - 13) 



-i 



deformability of ([7]) and (Q and does not admit Mohr-Coulomb and Tresca as limits, 
but results in a simple expression. The performance of the deviatoric shape of the yield 




3D 30 



60 60 




Figure 14: Deviatoric section (|53|) : effects related to the variation of /3 (a) and n (b). 



surface is analyzed in Fig. [TU where the solid lines correspond to the limit of convexity, 



25 



(3 = 1 and n = 3. The curves reported in Fig. [3] (a) are relative to the values of (3 = 0, 1, 2, 
whereas for Fig. [3] (b) n takes the values {1, 3, 5}. 

A limitation of the yield surface described by eqns (j7J) and ({53"]) is that the deviatoric 
section cannot be stretched until the Rankine limit. This can be easily emended assuming 
for g{6) our expression (J9]) or that proposed by Willam and Warnke (1975) (see also 
Menetrey and Willam, 1995) 

2(1 - e 2 ) cos 9 + (2e - 1) [4(1 - e 2 ) cos 2 6 + he 2 - 4e] a5 
ffW ~ 4(l-e 2 )cos 2 # + (2e-l) 2 ' ( } 

where e g]0.5, 1] is a material parameter, yielding in the limit e — > 0.5 the Rankine 
criterion and the von Mises criterion when e = 1. 

It is already known that the deviatoric section of the yield surface corresponding to 
eqn ( 1541 remains convex for any value of the parameter e ranging within the interval 
]0.5, 1], so that — from Proposition 1 — the function ([6]) equipped with the definition ( 154|) 
of the function g(8) is also convex. 

As a final example, we can employ function g(9) defined by the expression proposed 
by Gudheus (1973) and Argyris et al. (1974) 

2k 

9{e) = i + fc+a-^coss^ (55) 

where k e] 0.777, 1] is a material parameter. 

Otherwise, we can act on the meridian function. For instance, we can modify a 
Drucker-Prager criterion — which again fits in the framework described by eqn flH])— 
obtaining a non-circular deviatoric section described by eqn OH]) 



F(o) = — r (p + c) + q cos 



7T 1 

~ o cos_1 (7cos36>) 
o 3 



(56) 



where c is the yield strength under isotropic tension and T is a material parameter, or 
by eqn ( 1531 

F(tr) = -T(p + c) + q[l + (3(1 + cos 39)] 1/n , (57) 



or by the Gudheus/ Argyris condition (1551) 

F((t) = -T (p + c) + ^(1 + k + (1 - k) cos36) . (58) 

It may be noted that the yield criterion ( 158]) has been employed by Laroussi et al. (2002) 
to describe the behaviour of foams. In all the above cases, Proposition 1 ensures that 
for the range of parameters in which the Haigh-Westergaard representation of the yield 
surface is convex, the yield function is also convex. 
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4.3 A note on the behaviour of concrete and a generalization 
of Proposition 1 

In the modelling of concrete there is some experimental evidence that the deviatoric 
section starts close to the Rankine limit for low hydrostatic stress component and tends 
to approach a circle, when confinement increases. This effect has been described by 
Ottosen (1977) through a model which does not fit the general framework specified by 
eqn ([6]) and can be written in our notation in the form 

F(<T) = Aq 2 + B^- + C-p, (59) 

where A > 0, B > and C < are constants and g{9) is in the form (Q with (3 = 0. 
The criterion is therefore defined by four parameters. 

The above-expression (1591) of the yield function suggests the following generalization 
of Proposition 1: 

Proposition 2: Convexity of the yield function 

F{tr) = Aq 2 + B^- + f(p), (60) 

where A and B are positive constants, is equivalent to 

f >0, & g* + 2g> 2 -gg">Q, (61) 
which in turn is equivalent to the convexity of the surface 



(62) 



in the Haigh-Westergaard stress space. 

Proof. Let us begin assuming that (1611) holds true. In this condition Proposition 1 
ensures that f(p) and q/g{0) are convex functions of cr, so that (16"U|) results the sum of 
three convex functions and its convexity follows. Vice- versa, failure of convexity of f(p) 
immediately implies failure of convexity of (160]) since p is independent of 9 and q. Finally, 
let us assume that condition (161^0 is violated, for a certain value, say 9, of 9. The Hessian 
of 

Aq 2 + Bq/g(9), (63) 

as a function of two components of deviatoric stress Si and S 2 , is given by eqn. 
summed to a constant and positive definite matrix 



3A 



2 1 
1 2 



+ B 



27 (g 2 + 2g> 2 -gg") 
4 q 3 g 3 



5f 



■S1S2 



-S1S2 
SI 



(64) 
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Considering now the Haigh-Westergaard representation, it is easy to understand that we 
can keep 9 = 9 fixed and change Si and consequently 52 so that Si/S 2 remains constant. 
In this situation, while g and its derivatives remain fixed, the quantity 




(65) 



in matrix (I64p tends to +oo when Si tends to zero. Therefore, violation of ( IBTi) ? cannot 
be compensated by a constant term and function (160]) is not convex. 



Proposition 2 provides the conditions for convexity of the Ottosen criterion. Moreover, 
the same proposition allows us to generalize our yield function (E])-([9j) adding a q 2 term 
as in the Ottosen criterion. This leads immediately to 



where f{p) is given by eqn ([7]). 

5 Conclusions 

In the modelling of the inelastic behaviour of several materials, the knowledge of a 
smooth, convex yield surface approaching known-criteria and possessing an extreme shape 
variation to fit experimental results may be of undoubted utility. In particular, the fact 
that a yield function can continuously describe a transition between yield surfaces typical 
of different materials is of fundamental importance in modelling the de-cohesion due to 
damage of rock-like materials and the increase in cohesion during forming of powders. In 
the present paper, such a yield function has been proposed, which is shown to be capable 
of an accurate description of the behaviour of a broad class of materials including soils, 
concrete, rocks, powders, metallic foams, porous materials, and polymers. Moreover, 
in order to analyze convexity of our function, we have provided certain general results, 
holding for a broad class of yield conditions, which permit to infer convexity of the yield 
function from convexity of the yield surface. 
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□ 



F(<t) = Aq 2 



+ Bq cos (3 cos 1 (7 cos 3$) + f(p), 



(66) 
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APPENDIX A. Values of material parameters employed to fit experimental 

results in Sec. 13.31 

The values of the material parameters of the proposed yield function ([6])- (H employed 
to model the experimental data reported in Figs. [5Jl9]and Figs. [TTlfT2l are listed in Table [5] 
below. 



Table 5: Values of material parameters employed to fit experimental results in Figs. [5]|n] and Figs. lTTl[T21 
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Pc/fc = 3.5 x 10 7 


c/fc = 0.08 


1 x 10 6 
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APPENDIX B. Covexity of function (1551) 



We show that the deviatoric section described by eqn ( 153|) is convex, within the range of 
parameters reported in Tab. HI To this purpose, a substitution of eqn ( 1531) into condition 
d37D yields 

acos 2 (3#) + 6cos(3#) + c > 0, 0e[O,7r/3], (B.l) 
where the coefficients a, b and c are: 

a = /3 2 ( n 2 -9), (B.2) 
6 = pn(l + p)(9-2n), (B.3) 
c = n 2 (l + /3) 2 + 9/3 2 (l-n). (B.4) 

Since condition (IB. II) is bounded by a parabola, it suffices to analyze the position of its 
vertex to obtain the values of parameters reported in Tab. HI 
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